Pairing in Asymmetrical Fermi Systems with Intra- and Inter-Species Correlations 
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We consider inter- and intra-species pairing interactions in an asymmetrical Fermi system. Using 
equation of motion method, we obtain coupled mean-field equations for superfluid gap functions and 
population densities. We construct a phase diagram across BCS-BEC regimes. Inclusion of intra- 
species correlations result in stable polarized superfluid phase on BCS and BCS sides of unitarity 
at low polarizations. For larger polarizations, we find phase separations in BCS and BEC regimes. 
A superfluid phase exists for all polarizations deep in BEC regime. Our results should be apply 
broadly to ultra-cold fermions, nuclear and quark matter. 
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Pairing in two-species Fermi systems with unequal 
population is of great current interest and importance 
across a wide range of fields and systems. Examples are 
unequal density mixtures of fermionic cold atoms [3, 0] ; 
arbitrarily polarized liquid 3 He; superconductors in ex- 
ternal magnetic field [1] , in strong spin-exchange field 0, 
S @|, or with overlapping bands [3, Q; isospin asym- 
metric nuclear matter [9( and dense quark matter ex- 
hibiting color superconductivity [Io| . Unequal density 
cold fermions serve as prototypical systems, providing an 
unprecedented window into exploring superfluidity with 
tunable repulsive and attractive interactions. These are 
attained by sweeping across with s- or p-wave Feshbach 
resonances, thereby allowing the study of fermion ground 
states in both BCS and BEC regimes. 

Among the outstanding questions in asymmetrical 
Fermi systems is the nature of the ground state in the 
BCS and BEC regimes, and whether the BCS super- 
fluid state can sustain any finite imbalance between the 
species. Thus, it is important to arrive at a plausi- 
ble phase diagram as a function of pairing interaction 
strength and species imbalance. Two-species systems 
are conveniently characterized as two pseudo-spin sys- 
tems. It is believed that the BCS ground state in a finite 
magnetic field h, is robust against spin polarizations for 
h ~ A, (A being the superconducting gap); beyond this 
it becomes unstable to a normal state. For equal popu- 
lation cold atom systems, there is theoretical agreement 
with experiments that find superfluid states in both BCS 
and BEC regimes with a "smooth crossover" around the 
"unitarity limit" (diverging singlet scattering length a s ). 

For systems with population imbalance, various theo- 
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retical scenarios have been proposed 
field calculations H, 13[ find the superfluid state to be 
unstable to phase separation into superfluid and normal 
states or a mixed phase in the BCS regime; a super- 
fluid state stabilizes however deep in the BEC regime. 
Currently there is intense experimental efforts in unequal 
density cold fermion atoms. One experiment Q observed 
a transition from a polarized superfluid to phase separa- 
tion at a polarization ~ 10% near unitarity on the BEC 
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FIG. 1: (Color online) Polarization P vs s-wave coupling, 
— l/(kpa 3 ) phase diagram for asymmetrical fermions:(a) with 
a representative intra-species correlation strength g\ = 20, 
corresponding to l/(fcf.at) = 1.25; (b) without intra-species 
correlations, g\ — (also obtained in Ref [Hj]). Vertical line 
refers to the unitarity limit; PS1, PS2, PS to phase separated 
regions; N to normal state. 



side. To date, theoretical calculations have mostly con- 
sidered inter-species s-wave interaction, and have ignored 
intra-species correlations. Ho et al 1J] attempted to in- 
corporate triplet correlations i n a somewhat phcnomcno- 
logical manner; Huang et al 
implications for a FFLO state; 
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Monte Carlo calculations 
hint at a polarized superfluid phase near unitarity. 

In this paper, we address the issue of the nature of the 
zero temperature (T=0) ground state of an asymmetrical 
Fermi system for arbitrary repulsive/attractive interac- 
tion strength and polarization. We also examine if the 
BCS superfluid state can sustain a finite population im- 
balance. While the unequal density cold fermion systems 
may provide a way to test our results, our paper should 
have a broader appeal, viz. electronic superconductivity, 
nuclear and quark matter superfluidity, etc. Generally, 
both inter-species and intra-species correlations may be 
present in an asymmetrical Fermi system. These may 
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arise from the underlying fermionic potentials (atomic, 
electronic, nucleon-nucleon, quark-quark) or from effects 
of the medium, i.e. "induced" interactions [l?], EH- We 
include the simplest ones allowed by symmetry: s-wave 
contact interaction between the species, and a p-wave in- 
teraction within the species. Following Leggett and 
Eagles (20j , our discussion is in terms of BCS-type pair- 
ing in both the BCS and BEC regimes, with the chemical 
potentials for the two species determined self-consistently 
with the pairing gaps. We do a detailed stability analysis 
of the multitude of states obtained from our equations. 

Our findings are dramatically different from those with- 
out intra-species correlations, Our proposed phase dia- 
gram, Fig. la, shows that at T=0, for smaller polar- 
izations, and sufficiently large intra-species correlations, 
gapped polarized superfluidity (hereafter referred to as 
SFP) becomes stable on both BCS and BEC sides of the 
"unitarity" limit. Depending on the inter- species interac- 
tion strength, at some polarization, SFP becomes unsta- 
ble via a lst-order transition to phase separation, denoted 
by PS1. PS1 is characterized by a negative "susceptibil- 
ity", SP/Sh; P being the spin-polarization, and h the 
difference of the chemical potentials, playing the role of 
a "magnetic field". For a given intra-species interaction, 
and for sufficiently weak inter-species interaction, SFP 
and PS1 undergo transitions to the normal state on the 
BCS side. The gapped SFP persists into the BEC regime, 
sustaining progressively smaller polarizations. Deeper in 
the BEC regime, we find a superfluid phase (SF) at all 
polarizations. In the BEC regime, in addition to PS1, 
we find the existence of a somewhat different phase sep- 
arated state, PS2, characterized by positive "susceptibil- 



ity" , but not satisfying requisite superfluid ground state 
stability criteria. 

For our detailed study, we consider a two-species Fermi 
system with unequal "pseudo-spin" populations. To al- 
low for both inter-species and intra-species correlations, 
and noting that pseudo-spin rotation invariance would 
be broken by unequal chemical potentials, we adopt a 
pairing Hamiltonian given by 

Ho 
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where the pseudospin a =|, J, denote for example the 
two hyperfine states of ultracold Fermi atoms, cf is 
the fermion creation operator with kinetic energy ^ a = 
(-ko — Mo-! AV is the chemical potential of each of the 
species. g]} k ,, and gj^, are the interactions between the 
up and down spins respectively, and V is the volume. 
The singlet interaction, gjj,, is taken to be a constant, 
g . This is usually expressed in terms of s-wave scatter- 
ing length a s using (Ann 2 a s / m)^ 1 = g^ 1 + X) fc (2e*.) -1 . 
A mean-field decoupling is attained by introducing three 
order parameters or gap functions (a, a' =j, I) given by, 
A CTCT ' [k, q) = - YsH 9kk> < c -k'+ q /2aC k '+ q /2a' > This re- 
sults in a mean-field Hamiltonian given by: 
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We employ the equation of motion method using imag- 
inary time normal and anomalous Matsubara Green's 
functions, G crcr i(k, r), F aa i (k, r), respectively, and our 
mean- field Hamiltonian, H MF . The coupled equations 
in terms of A aa i are given by G (T(T '(fc,r) and F aa i{k,r): 

d T G aa \k,T) = -<S(r)<W -£fc+g/2a<w(fc,T) 

+ J2K»a(k,q)F (7 „ a ,(k,T) (3) 

d T F aa ,{k,T) = i- k +q/2aFaa{k 1 T) 

+ Y,K*»(k,Q)G*»cr'(k,T) (4) 



where r is the imaginary time variable. These equa- 
tions may be Fourier transformed in the usual way with 
t — » iw n , d T — > —iw n , where iw n = (2n + l)ir/j3 
are the Matsubara frequencies, n being an integer and 
/3 = l/k B T. 

Here we focus on a superfluid condensate of pairs with 
zero center- of -mass momentum, q. Thus, we do not con- 
sider the a ^ Fulde-Ferrel-Larkin-Ovchinnikov(FFLO) 
state [3, H|, but which may also be studied within this 
scheme. Solving the Fourier transformed equations at 
q = 0, we obtain the 2-point correlation functions: 

G„ A k,i Wn ) = W. + ^AA^-AA^) 
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where f a = A 2 b- a + A_ a b a - a- a b a b- a ; f aa = 
A 2 A- a +A tT A 2 _ a -A a a- tT b- a ; f a - a = A(A 2 + A (T A_ CT - 
a a b- a )(6 al - Sai); D = (A 2 + A T A X ) 2 + a^b i - 
^(A 2 ^ + A 2 a_ (T )6_ cr ; with a a = £k+q/2o - iw n , b a = 
-£-k+ q /2a - iw n - We have set A TT = A T ; A u = A x ; 
A = A. The excitation spectrum can be found by ex- 
amining the poles of the Green's functions, yielding the 
quasiparticle energies, 
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4A 2 [(f fcT - £ H ) 2 + (Ai - A 2 ) 2 ] . Various quanti- 
ties can now be obtained from our 2-point correlation 
functions. Thus, particle concentrations, n a for the two 
species (a =t, I) are given by 
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where A is even for I — +, and odd for I = — ; np(E k i) are 
the Fermi functions. Likewise the three gaps equations 
are given by (a, a' =f, j): 
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The above five equations are coupled, and can be solved 
self-consistently for the three gap functions, A,A|,Aj_ 
for either fixed particle concentrations, n^, nj,, or fixed 
chemical potentials, /if, fi^. 

We assume equal masses for the two species, and take 
the particle spectrum to be e k = h 2 k 2 /2m. We adopt 
standard definitions: polarization, P = (n\ — ri2)/(ni + 
n.2); mean chemical potential /1 = (fii + ni)j2; chemical 
potential difference h — (fii — /i 2 ) /2; Fermi momentum 
k Fa = (67r 2 ^)V 3 . Since £ fc /(*) -> ~> 

Jq 00 k F f(k/kp) d m&f^ , we can scale quantities having 
dimension of energy to e^- The inter-species interac- 
tion g G is expressed in terms of coupling constant r/ = 
— l/(kFa s ). For the intra-species triplet interaction, we 
take the separable form g%%, — giOj(k)uj(k')Y w (k)Yi (k') 1 
where we have taken <?ff = g^ = g~\. More generally the 
m = ±1 terms would also be present; however this choice 
allows us to explore the consequences of intra-species cor- 
relations while keeping the calculations tractable. As 
a check, we also consider different types of momentum 
dependence: (i) u>(k) oc const; (ii) ui(k) oc fcofc/(fc 2 + 
k 2 ), a generalization of the Nozieres and Schmitt-Rink 
scheme (HI; ( m ) u(k) oc exp [— (k/ko) 2 ], a Gaussian in- 
teraction; k a being a cut-off momentum; these give qual- 
itatively similar behavior. The first two forms of in- 
teraction require regularization due to ultraviolet diver- 
gence, while the third does not. With regularization, g\ 
can be expressed in terms of a triplet scattering volume 
a t 0: (^klat/m)- 1 = g- 1 +J2 k w(k) 2 /(2e k ). Thus, 
(3n/2ep)gi = g\ in our plots can be easily expressed in 
terms of at; e.g. g\ = 20 corresponds to l/(k F a t ) — 1.25. 



For arbitrary inter-species s-wave and and intra-species 
p-wave pairing interactions, and population imbalances, 
we obtain self-consistent solutions of the T = gap func- 
tions, A, Aj,AJ,, and chemical potentials, On the 
BCS side, for a given g a , the jj, gap A decreases with 
increasing intra-species interaction strength gi, while at 
the same time both Af,A^ (Aj 7^ AjJ increases, cross- 
ing at some value of g±. The suppression of A is more 
pronounced at larger polarizations. 

A proper construction of the asymmetrical Fermi sys- 
tem phase diagram requires a determination of stable 
ground states out of the manifold of paired condensates 
given by our equations ID, 23} . Accordingly, we carefully 
consider the stability criteria. The mean-field ground 
state energy as a function of the gaps at different po- 
larizations, P is given by: 



E G (A 1 A ] .A l ) = (m MF \V) 
= E o + Y / [^kaG aa {k,T^0- 

ka 
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where E Q = -|A T | 2 /.gi - |A x | 2 / 5 i - \A\ 2 /g D . To find the 
stability of the polarized superfluid state SFP, we con- 
struct the 3x3 stability matrix out of all partial deriva- 
tives = A, AjjAjJ, and check for posi- 
tive definitcness of the determinant of the matrix, and 
of all its upper-left sub-matrices. We supplement this 
with analysis of the "susceptibility" dP/dh. Thus, for 
a given 171, the stable polarized superfluid state SFP, in 
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FIG. 2: (Color online) Polarization P = ^| + "| vs 7l//_t = 

MT ; M | (/i acts like a "magnetic field") for different intra- 
species strengths gi at a fixed inter-species coupling r\ = 0.25. 



both BCS and BEC regimes, is characterized by Eq with 
a global minimum at non-zero gaps and self-consistently 
determined values of fj, a 's, and §77 > 0. SFP sustains 
larger polarizations for progressively larger g\ . 

Similar to the case without intra-species correlations, 
gx, there exists a maximum polarization, P max on the 
BCS side, beyond which we find no solution to the 
coupled equations; this determines the SFP/PS1 <-> N 
boundary (Fig. la). P m ax is slightly decreased at uni- 
tarity by g\. It decreases with increasing 77 = — l/kpa s . 
For a fixed gx, close to both BCS and BEC sides of uni- 
tarity, and for small polarizations, unlike-spin pairing has 
appreciable value in the polarized superfluid SFP. How- 
ever away from unitarity on the BCS side, At i decreases, 
and || and J, J, pairing becomes more dominant in SFP, 
as inter-species interaction becomes relatively weak com- 
pared to intra-species interaction. On the BEC side away 
from unitarity, on the other hand, Atj_ is more dominant, 
and with At and Aj, becoming negligible, SFP becomes 
unstable to phase separation, PS2. A SF phase emerges 
deep in the BCS regime with predominantly unlike-spin 
pairing at low polarizations and majority spin pairing at 
higher polarizations. 

The region PS1 in Fig la is characterized by negative 
"susceptibility" , ^£ and doe not satisfy the stability ma- 
trix criteria. For a given P and h, Eg is a maximum 
at the non-zero gap solutions, separating two local min- 
ima - a feature of phase separation into a normal and a 
superfluid component by lst-order phase transition. In 
this context, it is instructive to study P as a function 
of h/n for different values of gx, for a fixed 77 (shown in 
Fig.2 for 77 = 0.25 (BCS side)). For gx = gi « 15, the 
slope is vertical {gi is the value of (71 at which maximum 
polarization, P max occurs for a given rf). For g± < g\, 
the slope is negative, corresponding to the BCS super- 



fluid state being unstable to the normal state for h > A, 
but robust against polarization for h < A. For gi > g\, 
the singlet superfluid state can sustain a finite polariza- 
tion, which exhibits a behavior over the range given by: 
P oc ah + bh 3 + c; a,b,c being constants. The linear 
behavior is achieved for larger values of gx, and at low 
polarizations. In examining P vs h/fj, behavior beyond 
Pmax, we find, for a given gx, two solutions of P corre- 
sponding to one value of h//i. To make a connection to 
Fig. la, obtained for g\ — 20, we note that the allowed 
range of polarizations for SFP in P vs hj \x considerations 
corresponds to the 77 = 0.25 vertical line, terminating at 
the SFP-PS1 boundary (Fig la). The same line extended 
from SFP-PS1 boundary to PSTN boundary correspond 
to the polarization range bounded by the two solutions 
ofPat a/i/ju in Fig.2. The region PS2 in BEC regime, 
though characterized by ^ > 0, is not a stable super- 
fluid phase, since the stability matrix condition cannot 
be satisfied. The line separating PS1 and PS2 is proba- 
bly a metastable line, the position of which depends on 
the p-wave interaction strength. 



In summary, we find that the inclusion of intra-species 
correlations in asymmetrical Fermi systems results in a 
stable polarized superfluid phase SFP at low polariza- 
tions on both BCS and BEC sides of unitarity. We have 
discussed the nature of the paired states and transition 
to phase separated states. The SF phase obtained in the 
deep BEC regime in the case without intra-species corre- 
lations, also emerges here with dominant unlike-species 
pairing, accompanied by weaker majority-species pairing. 
Our results should be of broad interest as it should be of 
relevance to any asymmetrical Fermi system, with proper 
choice of interaction parameters. Here, our choice of pa- 
rameters appear to agree with cold atom experiment Q 
that found a SF to PS transition around ~ 10% polariza- 
tion around unitarity on the BEC side. Also, the max- 
imum polarization ~ 70% at unitarity is in agreement 
with experiments. Further experiments at low polariza- 
tions on both sides of unitarity are needed to test our 
detailed results. Experiments that measure differences 
in momentum distributions of two species could provide 
further test. Finally, our phase diagram indicates a tri- 
critical point (SFP,PS1,N phases) at low polarization on 
the BCS side, in addition to one on the BEC side at 
P ~ 1. This should lead to interesting study of the evo- 
lution of these two tricritical points at finite-T; we are 
exploring these effects. 



We would like to thank D. Allendcr, K. Bedell, J. En- 
gelbrecht, S. Gaudio, R. Hulet, and M. Widom for fruitful 
comments and discussions. We also acknowledge support 
from ICAM. 



■5 



[1] M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and 

W. Ketterle, Science 311, 492 (2006). 
[2] G. B. Partridge, W. Li, R. I. Kamar, Y. A. Liao, and 

R. G. Hulet, Science 311, 503 (2006). 
[3] G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963). 
[4] A. I. Larkin and Y. Ovchinnikov, Sov. Phys. JETP 20 

(1965). 

[5] P. Fulde and R. A. Ferrell, Phys. Rev. A 135, 550 (1964). 
[6] S. Takada and T. Izuyama, Prog. Theo. Phys. 41, 635 
(1969). 

[7] H. H. Suhl, B. T. Matthias, and L. R. Walker, Phys. Rev. 

Lett. 3, 552 (1959). 
[8] J. Kondo, Prog. Theo. Phys. 29, 1 (1963). 
[9] A. Sedrakian and U. Lombardo, Phys. Rev. Lett. 84, 

602 (2000); U. Lombardo, P. Nozieres, P. Schuck, H.- 

J. Schulze, and A. Sedrakian, Phys. Rev. C 64, 064314 

(2001). 

[10] M. Alford, J. A. Bowers, and K. Rajagopal, Phys. Rev. D 
63, 074016 (2001); J. A. Bowers and K. Rajagopal, Phys. 
Rev. D 66, 065002 (2002); T. Schafer and F. Wilczek, 
Phys. Rev. D 60, 074014 (1999). 

[11] P. F. Bedaque, H. Caldas, and G. Rupak, Phys. Rev. 
Lett. 91, 247002 (2003); H. Caldas, Phys. Rev. A 69, 
063602 (2004); M. M. Forbes, E. Gubankova, M. V. Liu, 
and F. Wilczek, Phys. Rev. Lett. 94, 017001 (2005); 
W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 
(2003); K. B. Gubbels, M. W. J. Romans, and H. T. C. 



Stoof (2006). 

[12] C.-H. Pao, S.-T. Wu, and S.-K. Yip, Phys. Rev. B 
73, 132506 (2006); Q. Chen, Y. He, C.-C. Chien, 
and K. Levin, cond-mat/0608454 (2006); Z.-C. Gu, 
G. Warner, and F. Zhou, cond-mat/0603091 (2006); 
M. M. Parish, F. M. Marchetti, A. Lamacraft, and B. D. 
Simons, cond-mat/0605744 (2006); D. E. Sheehy and 
L. Radzihovsky, Phys. Rev. Lett. 96, 060401 (2006). 

[13] D. E. Sheehy and L. Radzihovsky, cond-mat/0608172 
(2006); D. E. Sheehy and L. Radzihovsky, cond- 
mat/0607803 (2006). 

[14] T.-L. Ho and H. Zhai, cond-mat/0602568 (2006). 

[15] X. Huang, X. Hao, and P. Zhuang, cond-mat/0610610 
(2006). 

[16] J. Carlson and S. Reddy, Phys. Rev. Lett. 95, 060401 
(2005). 

[17] A. Bulgac, M. M. Forbes,and A. Schwenk, Phys. Rev. 

Lett. 97, 020402 (2006). 
[18] K.F. Quader,and K.S. Bedell, J. Low Temp. Phys. 58, 

89 (1985). 

[19] A. J. Leggett, J. Phys.(Paris) Colloq. 41 (1980). 
[20] D. M. Eagles, Phys. Rev. 186, 456 (1969). 
[21] P. Noziers and S. Schmitt-Rink, J. of Low Temp. Phys. 
59, 195 (1985). 

[22] M. Iskin and C.A.R. Sa de Melo, cond- 

mat/0602157(2006). 
[23] L. Viverit, C. J. Pethick, and H. Smith, Phys. Rev. A 

61, 053605 (2000). 



